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Abstract 

It is argued that in the case of a smooth transition across a (dilaton-driven) curva- 
ture bounce the growing mode of the vector fluctuations matches continuously with a 
decaying mode at later times. Analytical examples of this observation are given both in 
the presence and in the absence of fluid sources. In the case of multidimensional bounc- 
ing models the situation is different, since the system of differential equations describing 
the vector modes of the geometry has a richer structure. The ampliflcation of the vector 
modes of the geometry is specifically investigated in a regular five-dimensional bounc- 
ing curvature model where scale factors of the external and internal manifolds evolve 
at a dual rate. Vector fluctuations, in this case, can be copiously produced and are 
continuous across the bounce. The relevance of these results is critically illustrated. 



1 Formulation of the problem 



It is known that in conventional inflationary models vector fluctuations of the metric are 
not amplified by the pumping action of the gravitational field [1]. On the contrary both 
scalar and tensor modes of the geometry are copiously generated in the transition from a 
quasi-de Sitter stage of expansion to a radiation-dominated evolution. 

Vector modes of the geometry may have interesting phenomenological implications. For 
instance, as argued long ago by Wassermann [2] (see also [3]) the existence of large-scale 
vorticity (sometimes connected with a large-scale magnetic field) may have relevant implica- 
tions for the origin of global rotation of spiral galaxies. The possible existence of primordial 
vorticity has been the subject of various investigations. Harrison [4] (see also [5] for a 
complementary perspective) noticed that if vorticity is present over sufficiently large length 
scales in a plasma of ions and electrons, magnetic fields can be induced over scales compara- 
ble with the typical scale of the vorticity. The dynamical origin of primordial vorticity has 
been, since then, a subject of debate. In [6] it was suggested that cosmic strings with small 
scale structure may be responsible for the dynamical origin of vorticity. The dynamical 
friction between cosmic strings and matter may also be responsible for further sources of 
vorticity [7]. 

In spite of the fact that the considerations quoted so far may indeed produce large-scale 
vorticity, it is useful to take another approach and consider the possibility that, thanks 
to the pumping action of the gravitational field, the vector fluctuations of the geometry 
will be amplified, providing, at some later epoch, the initial condition for the evolution 
of the rotational fluctuations of the fluid. The motivation of the present investigation is 
to discuss the evolution of vector modes of the geometry in the case of a bouncing back- 
ground. As is known, the ten-dimensional degrees of freedom of a four-dimensional metric 
can be decomposed, with respect to three-dimensional rotations, into four scalar modes, 
two tensor modes (described by a divcrgencelcss and tracclcss rank-two symmetric tensor 
in three dimensions), and four vector degrees of freedom (described by two divergenceless 
vectors in three dimensions). Out of these ten degrees of freedom, four can be fixed by a 
choice of coordinate system. While the tensor modes are invariant under infinitesimal gauge 
transformations, the scalar and vector modes are not. 

While this work was in progress, Battefeld and Brandenberger [8] argued that it is not 
unlikely that vector modes of the geometry may be amplified in the case of four-dimensional 
curvature bounces as the ones provided by cyclic/ekpyrotic scenarios [9] or in the context 
of the pre-big bang scenario [10]. The authors of Ref. [8] worked in a four-dimensional 
conformally flat FRW model and discussed the dynamics in the context of an Einsteinian 
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theory of gravity. They found a growing-vector zero mode and correctly pointed out that 
the specific fate of such a mode is rather sensitive to the detailed dynamics of the bounce. 

There are, in principle, different bouncing dynamics that can be analysed. The simplest 
case is the one of a four-dimensional background without fluid sources but with a scalar 
degree of freedom (the dilaton). In this case, during the contraction the vector zero mode 
increases, but as will be argued in specific examples, its fate depends upon the dynamics of 
the curvature bounce and, in all the analytical examples reported in the present paper, it 
matches with a decaying mode after the bounce. A similar situation also arises, also in four 
dimensions but in the presence of fluid sources and together with a dilatonic background. 

In the multidimensional case the situation becomes increasingly interesting. There more 
scalar and vector modes arise and their quantum mechanical fluctuations can be amplified. 

For the sake of completeness, even if the scalar modes of the geometry are not the 
subject of the present investigation, we want to mention that curvature bounces have been 
the subject of various debates concerning the fate of the scalar fluctuations of the geometry 
(see [11] for a partial list of references). 

The plan of the present paper will be the following. In Section 2 the evolution of 
vector modes of the geometry will be analysed in the context of four-dimensional curvature 
bounces. In the string frame description, it will be shown that the vector modes present 
prior to the curvature bounce match continuously with a decreasing vector mode after 
the bounce. Specific analytical examples will be presented both in the absence and in 
the presence of fluid sources. In Section 3 the specific case of multidimensional curvature 
bounces will be discussed. The evolution of the vector fluctuations of a five-dimensional 
curvature bounce will be derived in fully gauge- invariant terms. In Section 4 the spectrum of 
the vector modes will be derived and it will be shown that vector modes can be amplified in 
the five-dimensional model of curvature bounce. Section 5 contains the concluding remarks 
with a critical summary of the main suggestions stemming from the present study. Finally, 
the appendix contains useful technical results needed for the derivation of the evolution 
equations of vector modes in more than four dimensions. 



2 Four- dimensional curvature bounces 



Consider, to begin with, the low-energy gravi-dilaton action in four space-time dimensions 
and in the string frame description [12]: 
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The potential V{ip) is, for the moment, a local function of the dilaton field ip and is the 
string length. A linear combination of the equations of motion derived from Eq. (2.1) leads 
to the well known form of the low-energy beta functions 



1 dV 



(2.2) 



In a conformally flat FRW background and assuming the dilaton field is homogeneous 

ds'^ = Gapdx'^dx^ = a'^{v)[dri'^ - dx% ip = ip{ri), (2.3) 
the component of Eq. (2.2) leads to the following relation 



2 dip 



(2.4) 



In Eq. (2.4) the prime denotes a derivation with respect to the conformal time coordinate 
r] while H = (In a)' and ip' = ip' — 37i is the time derivative of the so-called shifted dilaton, 
which is invariant under scale-factor duality. 

In this framework, the vector modes of the geometry are parametrized by two vectors, 
Qi and Wi, which are diver genceless, i.e. diQ^ = diW^ = ^: 



SGoi = -a'^{rj)Qi, 

6Gij = 2a^{r])d^iWj)=a\r])[diWj + djWi 



(2.5) 



Prom Eq. (2.5), the components Ricci tensor can be computed to first order in the fluctu- 
ations of the metric : 

6Ro^ = {W + 2n^)Q, - ^V^lQi + Wl], 

6Rij = -d^iQ'j) - ^w;') - 2{n' + 2re)d^iWj) - 2n[d^iW'^) + (2.6) 



Finally, using Eq. (2.6), the first-order fluctuation of Eq. (2.2) leads to the following 
pair of differential equations: 



H' -H'^ -nip' + 



a^dV_ 



Q^-lv^{Q^ + Wl) = 0, 



-5(,Q;.) - ^M^/) + {p' + H)%Q,) + ^l^j) 

n'-re-Wp' + --_ 



-2 



d(^,) = 0. 



(2.7) 



(2.8) 



'^In the present investigation the attention will be confined to spatially flat manifolds. However, our 

results can be extended to the case of non-flat spatial geometry by replacing, when needed, the ordinary 
derivatives with covariant derivatives with respect to the spatial part of the geometry. 
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Recalling the specific form of Eq. (2.4), it can be appreciated that various terms cancel in 
Eqs. (2.7) and (2.8). Defining now the new variable 



Vi = Qi + Wl, (2.9) 

Eqs. (2.7) and (2.8) can be written as 

VVi = 0, (2.10) 

V^ = iTp' + H)Vi. (2.11) 

The variable Vi is invariant under infinitesimal coordinate transformations, preserving 
the vector nature of the fluctuation. Indeed, under infinitesimal diffeomorphisms Qi and 
Wi transform as 



Wi = Wi + Ci, 

Qi = Qi- Ci, (2.12) 

where diC = is the gauge function. Hence, as anticipated, Vi is invariant under infinites- 
imal coordinate transformations. This observation also implies that the condition Wi = 
fixes the gauge freedom completely and is analogous, in this sense, to the conformally New- 
tonian gauge often employed in the analysis of scalar perturbations ^. This gauge choice 
would have led to the evolution equations for Qi, coinciding, in such a gauge, with the 
evolution of the gauge-invariant vector Vi. We will go back to the gauge issue during the 
discussion of vector fluctuations in the multidimensional case. 

Prom Eqs. (2.10) and (2.11) we can deduce that only the zero mode is dynamical. In 
this case its evolution can be written as 

Vi = a a{ri) (2.13) 

where Cj is an integration constant. During a conventional pre-big bang phase [13], as 
ry — >^ 0~, the consistent solution of Eq. (2.4) and of the other equations of the background 
can be written, in the absence of dilaton potential, as [13] 

1 



^Notice that, in the scalar case, the longitudinal gauge condition (imposed on the gauge parameters 
preserving the scalar nature of the fluctuation) eliminates the off-diagonal entries of the perturbed metric. 
In the case of vector modes the situation is, in a sense, opposite since the gauge choice fixing the coordinate 
system completely keeps only the off-diagonal elements of the perturbed metric tensor. 
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where t is the cosmic time coordinate obeying the usual differential relation a{r])dri = dt. 
Inserting Eq. (2.14) into Eq. (2.13) 



V3+1 

\t\ V-s . 



(2.15) 



This growing mode diverges in the limit t — 0~. This is essentially one of the arguments 
put forward in Ref. [8], with the difference that, in the present case, the discussion has 
been conducted, from the very beginning, in the string frame description. 

The aim will now be to show that in a simple model of low-curvature transition, the 
growing vector mode will match to a decreasing solution in the post-bounce epoch. As was 
shown in [10] a smooth transition at low curvatures can be achieved if the dilaton potential 
is a local function of Tp. The generally covariant action describing this model is given, in D 
space-time dimensions by 



(2.16) 



S = --^ / d^'-^x^le-^ \r + G'^^Va^Vpip + V{ip) 

where Sm represents the contribution of effective (fluid) sources and where 

V = Vie-^), 

e"^^"^ = / d''^\/\GM\e-''^''^^G-f^dM'^)dM'^)^{^i^) ' </'(^)) -(2.17) 

The variation of the action (2.16) with respect to and (p leads to the following equations: 



1 



G'^'^Va'pVpip - 2G''f^VaVpp - V 



R + 2G'^NoNpP> - G'^'^da'pdpip + V 



dV 
dtp 



where ^ 



and 



d^(pdy<p 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



^In Eq. (2.19) as well £is in Eqs. (2.22) and (2.23) the prime denotes a derivation with respect to the 
argument of the given functional and not the derivative with respect to the conformal time coordinate. The 
two notations cannot be confused since the latter is only employed when dealing with the explicit form of 
the equations on a given background. 
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is the induced metric. In Eqs. (2.18) and (2.19) the following integrals 



^ / d''w^\G{w)\^G-PdMw)dpip{w)5\ip{x) - ^{w)), (2.22) 



have also been defined. 

By combining Eqs. (2.18) and (2.19) in order to eliminate the Ricci scalar we obtain 
the following equation 



dV 
dip 



G 



111/ 



2^1 — ^T)iti/. 



(2.23) 



If Ji and I2 are set to zero and if dV/dip is replaced by dV/dip, Eq. (2.2) can be formally 
recovered from Eq. (2.23). 

In the case of a homogeneous dilaton and for a conformally fiat metric of FRW type 
Eqs. (2.18) and (2.23) lead, in D = 4 and in units 2A5 = 1, to the following system of 
equations: 



Tp^ -?,H'^ -V = e'^p, 

H = TpH+ ^e^p, 

.. .2 9 dV 
2lp-Tp -dH^ + V - — = 0, 

dip 



(2.24) 
(2.25) 

(2.26) 

(2.27) 



(2.28) 



p + dHp = 0, 
where the stress tensor of the fluid sources 

t; = (p+p)ux-p'^;: 

has been assumed and where the overdot denotes a derivation with respect to the cosmic 
time coordinate. In Eqs. (2.24)-(2.26) the energy and pressure densities have been rescaled 
as 

p = a^p, p = a^p. (2.29) 
In the case p = p = 0, Eqs. (2.24)-(2.27) admit the solution [10] 

V{Tp) = -Foe^^, (2.30) 

(2.31) 



a{t) = ao 



T + \/r2 + 1 



1/V3 



(2.32) 
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where 

t e~'^'^° 

The scale to fixes the width of the curvature bounce ^. 

Another solution, which is effective in illustrating some aspects of the present analysis, 
may be obtained if the scale factor and the energy density of the fluid are both constant. 
In this toy example, Eqs. (2.24)-(2.27) may be integrated easily and the result is 

a = l, p = po, P = 0, 
V = -Vbe^'^, e'^ 



i + (W 



e 



4 



Po = T2 =1^oe"^°- (2-34) 



The evolution equations of the vector modes can be rederived by perturbing to first 
order Eq. (2.18). In the case without sources, the result coincides in = 4, with Eqs. 
(2.10) and (2.11). In the case with sources Eqs. (2.10) and (2.11) become 

Wi = -a2e^(p + p)Vi, (2.35) 
V( = {T^ + H)Vi, (2.36) 

where Vj is the fluctuation of the velocity field and where we used the explicit perturbed 
form of the energy— momentum tensor 

6Tf = {p + p)Vi, «° = ^' 5ui = aVi. (2.37) 

Equations (2.30)-(2.32) can be used in order to obtain the explicit solution for the 
gauge-invariant vector mode, which is 

V,^,Ml±^^f^, (2.38) 

where, as before, c,- is an integration constant carrying vector indices. This shows that for 
t <C —to, Vi ^ t , as obtained previously. However, for t ^ to the dangerous growing 
mode turns into a decaying mode, i.e. ^ 

f{i-V3)/V3_ Furthermore, no divergence is 

present for t ^ 0. 



*The possibility of expressing the evolution of the background geometry in terms of the rescaled cosmic 
time T will also prove useful in the numerical discussion of the spectrum of vector fluctuations in multidi- 
mensional models. 
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Figure 1: The evolution of the vector zero mode in the case without fluid sources (full line) 
and with fluid sources (dashed line). 

In the case of the solution given in Eq. (2.34), Eqs. (2.35) and (2.36) lead to 

Vi = a- 



'l + r2' 



Vi = ^, (2.39) 
oopo 

where k denotes the Fourier mode. The time evolution of the gauge-invariant vector mode 
is illustrated in Fig. 1 for the two results obtained in Eqs. (2.38) and (2.39). 

The analysis reported so far shows that in some simple models of curvature bounce the 
growing mode in the vector sector matches with a decaying mode after the bounce. 



3 Multidimensional curvature bounces 



The analysis of gauge-invariant fluctuations in multidimensional and/or anisotropic cosmo- 
logical models has been discussed, through the years, in different contexts. One of the first 
attempts to adapt the original Bardeen [14] formalism to a multidimensional framework 
was discussed in [15] with the aim of performing a general discussion of the main features of 
cosmological fluctuations in Kaluza-Klein models [16]. Of related interests are some works 
discussing the theory of cosmological fluctuations in Bianchi models [17, 18]. In [19], a 
perturbative framework for the gauge-invariant discussion of multidimensional pre-big bang 
models was introduced. 
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The problem of gauge-invariant fluctuations in models with extradimensions received 
new momentum in the context of brane models, and different approaches (see, for instance 
[20, 21, 22] and references therein) were developed in order to study the possible localization 
of modes of various spin of both gravitational and non-gravitational origin. In [23] vector 
modes of the geometry were analysed in the case of a six-dimensional Abelian vortex leading 
to the localization of gravitational interaction (see also [24] for a more general discussion of 
codimension-two braneworlds and of their vector excitations). 

It is useful to notice, preliminarily, that the action given in Eq. (2.16) also admits 
classical solutions in the presence of internal (contracting) dimensions. In this section, 
the vector fluctuations of multidimensional cosmological models will be studied. The line 
element will be written, in this CclSGj clS 

ds'^ = Gf^^^dx^dx" = [dt'^ - a^{t)gijdx'dx^] - h^{t)gABdy^dy^ , (3.1) 

where G^jj/ is the D-dimcnsional metric; gij{x) and gAsiy) arc the metric tensors of two 
maximally symmetric Euclidian manifolds. The indices yU, v run over the full D = d + 
n + 1-dimensional space-time; the d-expanding dimensions are typically three and the n- 
contracting dimensions (labelled by capital Latin letters) will be, for a while, taken to 
be generic but the case n = 1 will be of particular interest in the following part of the 
investigation. 

From Eq. (3.1) in the particular case gij — 6ij and gAB — ^AB^ Eqs. (2.18) and (2.19) 
can be written in explicit terms by introducing the expansion rates H = a/a and F = b/b 
(the dot denotes a derivation with respect to the cosmic time coordinate): 



^ -dH"^ -nF^ -V = (3.2) 

H = Hip, (3.3) 

F = FTP, (3.4) 

Ot r 

2ip-i^-dH^-nF^ + V- — = 0. (3.5) 

op 



Equations (3.2)-(3.4) come, respectively, from the (0,0), {a,b) components of Eq. 

(2.23), while Eq. (3.5) can be derived from Eq. (2.19). 

In the case of V = — Vbe^'^, Eqs. (3.2)-(3.5) are solved by the following ansatz 

a{t) = (r-F + l)vfc, b{t) = {r + Vr^ + (3.6) 
p = po-\ln[l + T% r = ^. (3.7) 

The solution given in Eqs. (3.6) and (3.7) describes the expansion of the d dimensions and 
the contraction (at a dual rate) of the n internal dimensions. 
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Recalling the relation between the cosmic time coordinate t and the conformal time 
coordinate rj, i.e. a{r])dr] = dt, Eqs. (3.2)-(3.5) can be written as 

Tp'^ - dn"^ - nT^ - Vo? = 0, (3.8) 
n' = 'H^ + 'H^, (3.9) 
= nr + TT^, (3.10) 

r\T r 

^2 _ + + dn^ + nj^2 + — - Fa^ = 0. (3.11) 

cRf 

where, in full analogy with the four-dimensional case, the rates of expansion in conformal 
time have been defined: 

H = ilnay, J^=(ln6)', ' = (3.12) 

A particularly interesting case is the one where only one extradimension is present (i.e. 
n = 1). In this case, setting n = 1 and d = 3 in Eq. (3.6) the scale factors are 

«(0 = ^t)= Vr + Vr^ + l. (3.13) 

Prom Eq. (3.13), it can be easily deduced that the scale factor a{t) smoothly interpolates 
between a{t) ~ t~^^^ for t — >■ — oo and a{t) t^^"^ for t +oo. Moreover, b{t) will contract 
like for t <ti —tQ and will contract like t^^/^ for t ^ to. Both the scale factors and the 
curvature invariants are regular in the origin and for \t\ oo. Finally, a pleasant feature 
of the case d = 3 and n = 1 is that in the limit t +oo the usual radiation-dominated 
evolution is correctly reproduced. In Fig. 2 the time evolution of the scale factors a(r) and 
6(r) and of the associated rates, i.e. H{t) and F{r), is illustrated. This model describes 
the decoupling of internal and external dimensions. 

In the following part of the present Section, the vector modes arising in the model of 
Eq. (3.13) will be specifically analysed. The equations of motion will be derived in a gauge- 
independent way. The fluctuations of a five-dimensional line element of the type of Eq. 
(3.1) are, in general, described by 15 degrees of freedom which can be formally decomposed 
in terms of their transformation properties with respect to three-dimensional rotations as 



SGoo 


= 2a2(/., 




SGoi 


= -a'^diP 


- a^Qi, 


SGoy 


= ~abC, 




SGij 


= 2a^ 11)51 j - 


- 2a^didjE + a^hij + a^{diWj + djWi 


SGiy 


= —abdiD 


— abHi, 


SGyy 


= 2b% 





(3.14) 
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Figure 2: The time-evolution of the scale factors a(r) (solid line) and &(r) (dashed line) 
introduced in Eq. (3.13). The time evolution of the associated rates of expansion and 
contraction is also reported. 

where 

djh^i =0, Wi = 0, (3.15) 
diW' = 0, diQ' = 0, diW = 0. (3.16) 

According to (3.15) hij carries two degrees of freedom since it is symmetric in the two 
indices. Equation (3.16) tells us that Wj, Qi and Hi carry overall six degrees of freedom. 
The remaining degrees of freedom are seven scalars, which will not be specifically discussed 
since the attention of the present analysis is mainly centred on the vector modes. 

The vectors Wj, Qi and Hi obey a set of differential equations, which we ought to derive 
now. Notice, preliminarily, that in terms of the divergenceless gauge parameter C,i the 
vectors Wj, Qi and Hi transform, for infinitesimal diffeomorphisms preserving the vector 
nature of the fluctuation, as 

Qi — Qi ~ Ci: 

Wi = Wi + Ci, 

Hi = Hi- ^dyQ, (3.17) 

As in the four-dimensional case, gauge-invariant fluctuations can be defined, namely 

Vi = Q, + Wl, (3.18) 
Z, = Hi + ^dyWi. (3.19) 
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Notice that (3.18) corresponds to the gauge-invariant fluctuation of the four-dimensional 
case. 

In order to derive the coupled evolution of Vi and Zi, Eq. (2.18) shall be linearized to 
first order in Qi, Wi and Hi. By inserting the definitions (3.18) and (3.19) into the obtained 
equations, the system will only contain the relevant gauge-invariant combinations, i.e. Vi 
and Zi. 

From Eq. (2.18), the general form of the various components of the perturbed equation 
can be written, in a compact form, as 



o-^SQl - ^[W + d^Qi] - ^[diW^' + dmi] = 0, 



2 62 



a2 dV 



(3.20) 
(3.21) 
(3.22) 



In order to derive Eqs. (3.20)-(3.22) the perturbed connections reported in the appendix 
(see Eq. (A.3)) have been inserted in the linearized form of Eq. (2.18), bearing in mind 
that the fluctuation of the induced metric (2.21) implies 



'^70 = 



(3.23) 



Prom Eqs. (3.20)-(3.22), using Eqs. (A.7)-(A.9) reported in the appendix, the following 
expressions are obtained: 
1 



1 



{diQ^ + &iQi)' - (^ + 'K)(diQ^ + &>Q^ 

{diW^ + d^Wi)" + (^ + n){diW^ + d^i)' 



diW +&>Hi 



262^^ 



diW^ + dmi 



0, 



Qi 



dV 

Tp'^ + SH^ + - 2ip" + 270 - Va" + —a' 

cKp 



+ n)Hi 



y 
b_ 
2a 



+ i^dyV^Wi 



Q'i + {n-2J^-ip')Qi 
Hi' - (^ + n)H'i - V^Hi -{n- J'fHi - V^Hi 



0. 



(3.24) 



(3.25) 



(3.26) 



Recalling the definitions of the gauge- invariant fluctuations given in Eqs. (3.18) and (3.19), 
Eq. (3.24) becomes 



-±dy[d^Z^+d^Z,] + \ 



{diV^ + d^Vi)' - (^ + n){diV^ + d^Vi] 



0, 



(3.27) 
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By virtue of Eq. (3.11), the first term of Eq. (3.25) vanishes. The use of the gauge- 
invariant fluctuations Vi and Zi into Eq. (3.25) then leads to 

-^M^i + ^^IVi + ^VV, = 0. (3.28) 

Finally, inserting Eqs. (3.18) and (3.19) into Eq. (3.26), the last gauge-invariant evolution 
equation will be: the (iy) equation 

^dy[v^ + {H-2j^- ip')Vi\ - ^[z'Z - {ip' + n)z[ -{n- j^fz^ - v'^z,] = o. (3.29) 

Equations (3.27)-(3.29) can also be derived with gauge-dependent methods, as can be 
expected from general considerations. For instance, fixing the gauge in such away that Wi = 
the set of equations (3.27)-(3.29) can be exactly reproduced, as was checked by explicit 
calculation. By formally taking the limit Zj ^ and by setting to zero the derivatives with 
respect to the internal coordinate y, the four-dimensional equations derived in Eqs. (2.10) 
and (2.11) are formally recovered. 



4 Spectrum of vector modes 

For the purposes of this investigation it is interesting to consider the spectrum of the zero 
modes with respect to the internal dimensions. This is achieved by setting to zero all the 
terms containing one (or more ) derivatives with respect to y. Thus, from Eqs. (3.28)-(3.29) 
the following system: 

V^Qi = 0, (4.1) 
Q'i = {n + ip')Qi, (4.2) 
H'/ -{n + i^)H[ -(n- J^fHi - V'^Hi = 0, (4.3) 

is obtained. The evolution of Hi and Qi is, in this approximation, decoupled. Furthermore, 
it can be noticed that the evolution of Qi is the one already discussed in the four-dimensional 
case, with the minor difference that, in the present case, the background evolution differs 
from the four-dimensional case. Prom Eq. (4.2) we then have 



lr{T^ + 1) + (r2 + 1)3/2 
Qr{t) = -^^-1 — , (4.4) 

which also implies that Qi ~ t~^/'^ for t S> to- Hence, even if the time dependence is different 
from the four-dimensional case, it is still true that the Qi decreases at large (positive) times 
even if the solution is growing just before the bounce. 
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Let us now concentrate on the evolution of Hi whose equation can be written, in Fourier 
space and for a single vector polarization as^ 



u" + [k'^ - U{ri)]u = 0, 



-lp/2 \ // 



W2 



(4.5) 
(4.6) 



where Hi = ^Jae^l'^Ui. This equation will be numerically solved in a moment for different 
values of k. However, in order to cross-check the numerical results it is useful to anticipate 
the analytical estimate. For this purpose a useful observation is that, for times that are 
sufficiently far from the bounce, the asymptotic solutions are an excellent approximation 
for the background evolution. 

With this logic in mind, for t <C — io the solution given by Eqs. (3.6) and (3.7) can be 
approximated, in conformal time, by 



a(r?)~(-r7) ft(^) ^ (_^)i/3^ ^^__ln(-^). 

For t 3> Iq, the solution (3.6) and (3.7) behaves, in conformal time, as 

1 



a(r/) ~ r?, 6(r/) ~ - 



ip ~ — 21nr7. 

Therefore, in the two opposite limits, Eq. (4.5) is given, respectively, by 

41.2 _ y 



Ui + 



k' 



k' 



4^2 

4^2 _i 



4^2 



Uj = 0, 77 < -r/i, 
ttj = 0, ?7 > ?7i. 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



where r/i ~ r/o is a time scale comparable with the size of the bounce. 



Inserting Eqs. (4.7) and (4.8) into Eq. (4.6) the indices /x and v are determined to be 
= 2/3 and /x = 2. Equations (4.9) and (4.10) can be solved exactly in the two limits. 
Now, imposing quantum-mechanical initial conditions for 77 —00, 



U-{j]) 



e"2 



Hi'\-kv), 



■n < -m, 



(4.11) 



the solution for ij — -|-oo can be written in terms of the mixing coefficients c±{k) 



u+{ri) 



= v![!?gif (,.+1/2) 



c^{k)Hj^\kr)) + c+{k)e-'<i'+^/^^H(^\kri) 



ri>riQ, (4.12) 



^Notice that in the following the index denoting the Fourier mode will be omitted in order to avoid 
confusion with the vector polarization. 



14 



where H^J\z) = H(%{z) = J„{z) + iY^{z) are the usual Hankel functions of argument 
z and index a [25]. Notice that the phases in Eqs. (4.11) and (4.12) have been carefully 
selected in such a way as to recover mode functions with canonical quantum mechanical 
normalization in the limits ry — ^ ±00. 

The two solutions (4.11) and (4.12) can be continuously matched across the bounce, i.e. 
imposing «_(— r/i) = it+(r/i) and similarly for the first derivative. The mixing coefficients 
turn out to be 



c+(A:) = 4a;ie-('^+'^+i)/2 - + ^ + ^ gW(xi)gW(a;i: 



4 



xi 



^ + ^ + l//f(x,)//«(xi) 



C-{k) = i—xie 



Xl 



+<)(xi)if«i(xi) + iJ^i(xi)iJ«(xi) 



(4.13) 



(4.14) 



satisfying the exact Wronskian normalization condition |c+(/c)p— |c_(fc)p = 1. In the small 
argument limit, i.e. krji ~ hrjo <§C 1 the leading term in Eq. (4.14) leads to 

c-{k) ~ e^^(^-'^)/2^r''"''('^ + - l)r(/^)r(z/), (4.15) 

which implies, after inserting the specific values of /i and u that the number of produced 
vector quanta is given, for each polarization, by 

\c-{k)\^ ^ \k/ki\-^^, 7 = (^ + i^) ~ 2.6, (4.16) 

implying, in turn, that the energy density per logarithmic frequency interval, i.e. k'^\c-{k)\'^ , 
increases at large distance-scales (small k) as k~^'^'^. 

The same estimate will now be obtained through a totally different procedure, where 
the evolution equation for Hi is solved directly using the exact background solutions. Since 
the background solutions have a simple form in cosmic time, it is more practical to solve the 
evolution equation of the fluctuations directly in cosmic time. This will not affect the deter- 
mination of the mixing coefficients, which are independent of the specific parametrization 
of the time coordinate. 

Equation (4.3) becomes, in cosmic time, 

Hi - TpHi + [u^ -{H- Ff]Hi = 0, (4.17) 

where ojit) = k/a{t) is the physical momentum. By defining now the rescaled variable 
Vi = e-^/'^Hi, Eq. (4.17) becomes: 



0. (4.18) 
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Inserting Eqs. (3.6) and (3.7) into Eq. (4.18), the following simple equation for a single 
vector polarization can be obtained: 



+ 



3(t2 + 2) 



v = 0, 



(4.19) 



.tV^MH: 4(r2 + 1)2_ 
where t = t/to is the usual rescaled cosmic time and k = kto. 

Notice that the variable v defined in Eq. (4.19) is simply related to the variable u defined 
in Eq. (4.5): 

Hi = ^K^l'^u, = (Pl'^vi. (4.20) 
Since v = ^JHu, the quantum-mechanical initial conditions will be, in terms of v 



i J ujdt 



(4.21) 



for t ^ — oo. For t +oo the solution can again be written in terms of the mixing 
coefficients c± as 



v{t) = 



1 



(4.22) 



Equation (4.19) was solved for different values of k (or k = kto, the dimensionlcss comoving 
momentum). Quantum-mechanical initial conditions are given when the given mode is 
inside the horizon at the initial time of integration. From the numerical solution for the 
real and imaginary parts of v and v, i.e. 



v = f + iq, v = g + ip, 



(4.23) 



it is then easy to obtain the mixing coefficients after the relevant mode re-entered during 
radiation, i.e. 



|c+(A;)|^ + |c_(A:)|^ = - 



H 



H 



\c+ik)\^-\c.ik)\^ = 2{qg-fp). 



(4.24) 
(4.25) 



The results of this calculation are reported in Fig. 3. Recall that, according to Eq. (4.25), 
|c+(A;)p — |c_(A;)p = 1 since |c+(A;)p — |c_(A;)p is proportional to the Wronskian of the 
solution. From Fig. 3 it is apparent that for an increase of one order of magnitude in k, 
|c_(fc)p increases by 5.1 orders of magnitude, exactly as predicted in Eq. (4.16). 



5 Concluding remarks 

In the present paper the evolution of vector modes of the geometry has been discussed in the 
case of bouncing cosmologies both in four and in higher dimensions. The main suggestions 
following from a series of detailed examples can be summarized as follows: 
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Figure 3: The numerical calculation of the mixing coefficients for different values of k = ktQ, 
the comoving momentum rescaled through the size of the bounce. The dashed line is Eq. 
(4.25) whose right-hand side has been numerically computed. 

• in the four-dimensional case, as also noted in [8], a growing vector mode is present in 
the regime preceding the curvature bounce; 

• however, in the simple analytical models described in the present paper, the growing 
mode matches with a decaying mode after the curvature bounce, so that no further 
dangerous divergences are introduced; 

• in higher dimensions, more vector modes are expected and their spectrum was com- 
puted in a simple five-dimensional model of curvature bounce, where the compactifi- 
cation of the internal dimension is achieved dynamically; 

• the typical spectra of vector modes obtained in the specific (but fully consistent) 
example are red after the modes reenter the horizon. 

The most problematic aspect of this analysis is represented by the possible occurrence of 
red spectra, which may be rather dangerous in models based on curvature bounces. 
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A Vector modes of the geometry 



The three divergenceless vectors appearing in five dimensions are defined as 

6Gij =a\diWj+djWi), 

SGiy = -abHi (A.l) 



and 



91 

a? 



5G'y = ^. (A.2) 
ab 

Notice that Eqs. (A.l) and (A.2) are general in the sense that no specific gauge choice has 
been imposed. Prom Eqs. (A.l) and (A.2) 

sr^ = -^{diQj + djQi) - n{diWj + djWi) - ^{diW; + djWi), 

5T{^ = \{d,Qi - - \{&Wl + d,W^'), 

b ■ b'^ ■ 

= -l^yQi + ^[^^ + + •^)^^]' 
STl = ^^[H', + {J^-H)m-^dyQ, 



2 

= -nQ%j + ^d^diWj + djWi) - \dj{d^Wi + diW^) - \di{d^w^ + d^w^), 

ST\y = ^{diH^ - d^Hi) - \dy{d^Wi + diW^). (A.3) 

Prom the expressions of the Christoffel connections, the first-order fluctuations of the Ricci 
tensors can be determined: 
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+^dy[Hl + {J'-nm, (A.4) 
SRij = ^dyidiHj + djHi] - ^i^idiQj + djQiY + {2H + J'MQj + 5,Qo} 

-^id,w;' + djwn - \{2n + rmw'^ + d^w[) 

„ , 1 



+^dy{diWj + djWi) - -{diWj + djWi){2n' + + 2hj^), (a.5) 



5Riy = ~dy[Q[ + {AH - T)Qi\ - ^dyV^Wi 



+^[Hi + i^n + J^)H[ + {H' +J^' + H^ + 5HJ^)Hi - V^i/ij. (A.6) 

In Eqs. (A.4)-(A.6), denotes the Laplacian of the internal dimensions. The perturbed 
Ricci tensors with mixed indices are also useful for a swifter calculation of the perturbed 
equations: 



5RI = ^[2H' 2H^ - 2nT + + i^^yQ" + 

+ d^Wi)" + {2n + T){diW^ + d^Wi)'^ - ^dy[diW^ + d^w,], (a.s) 

^R'i = + m - ^)Qi] + ^9y^^Wi 

~^\^H'/ + {2H + J^)Hl + {r' -H' - 3H^ + 3HT)Hi - V^i^ij. (A.9) 

Clearly, since vector fluctuations are parametrized by divergence-less vectors in three dimen- 
sions, their contribution to the first order fluctuation of the Ricci scalar vanishes. Hence, 
recalling the definition of the Einstein tensor Q^, 

6gi = 6Rl sgl = 5Ri, sgf = sRl (a.io) 
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